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Abstract. Collisions between 248 Cm and 48 Ca are systematically calculated by time- 
dependent density functional calculations with evaporation prescription. Depending on the 
incident energy and impact parameter, fusion, fusion- fission, and quasi-flssion events are 
expected to appear. In this paper, the evaporation prescription is introduced, which is expected 
to be rather important to heavy-ion reactions producing superheavy nuclei, where the heavier 
total mass can be related to the higher total excitation energy. 



1. Introduction 

The synthesis of superheavy chemical elements [TJ [2] in the laboratory is one of the biggest 
challenges in nuclear physics. It is an attempt for clarifying the existence limits of all the 
chemical elements, as well as the completion attempt of the periodic table of chemical elements. 
We are concerned with heavy-ion collisions 

248 Cm + 48 Ca 

with different energies and impact parameters in this paper (Fig. [T]). Let A and Z be the mass 
number and the proton number, respectively. The neutron number N is defined by N = A — Z, 
so that N/Z of 248 Cm and 48 Ca are 1.58 and 1.40, respectively. If fusion appears, 296 Lv (= 
296 116) with N/Z = 1.55 is produced. 

Fast charge equilibration [3] is expected to appear in low-energy heavy-ion reactions with an 
incident energy of a few MeV per nucleon. It provides a very strict limitation for the synthesis of 
superheavy elements. Actually, the N/Z of final product is not above nor below the N/Z of the 
projectile and the target (1.40 < N/Z < 1.58 in this case) in the case of charge equilibration, so 
that the proton-richness of the final product follows. Although the actual value of N/Z depends 
on the two colliding ions, its value for the merged nucleus tends to be rather proton-rich for 
a given proton number of the merged system. This feature is qualitatively understood by the 
discrepancy between the /3-stability line and the N = Z-line for heavier cases. 

In this paper, based on the time-dependent density functional theory (TDDFT), the 
reaction dynamics of 248 Cm + 48 Ca is discussed focusing on both, thermal property and charge 
equilibration dynamics. Starting with specific features in TDDFT, an implementation of an 
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Figure 1. (Colour online) Potential energy as function of nucleus- nucleus distance R is shown 
for 248 Cm+ 48 Ca. The binding energies of 248 Cm, 48 Caand 296 Lv are calculated by static density 
functional calculations, where the former two nuclei are actually used in TDDFT calculations. 
The Coulomb potential is shown by the red curve. 

evaporation prescription in TDDFT is explained (Sec. [2]); the TDDFT calculations are presented 
for both, before and after the evaporation prescription in Sec. [3j and the reaction dynamics is 
summarized in Sec. [U 

2. Identification of evaporation residues 

2.1. Treatment of the thermal property 

Self-consistent time-dependent density functional calculations are employed in this paper. 
TDDFT reproduces the quantum transportation due to the collective dynamics. In this sense, 
what is calculated by the TDDFT can be regarded as products after several 10 -21 s, which 
corresponds to a typical time-scale of low-energy heavy-ion reactions (1000 fm/c), as well as 
to the inclusive time interval of any collective oscillations such as giant dipole resonance, giant 
quadrupole resonance and so on. Meanwhile, thermal properties such as the thermal instability 
are not directly taken into account in TDDFT. Indeed, the Skyrme type interaction used in 
TDDFT (for example, see Ref. [4j) is determined only from several densities. It is important 
that the most effective cooling effect arises from the emission of particles, and therefore it is 
expected that the break-up or fission of fragments including rather high internal excitation 
energy is suppressed in the TDDFT final products. The additional thermal effects leading to 
the break-ups of fragments should be introduced. 

Here is a fact that simplifies the treatment of thermal effects, that is, the difference of the 
time-scales. Different from the typical time-scale of low-energy heavy-ion reactions, the typical 
time-scale of the thermal effects is estimated by the typical time interval of collision-fission 
(fission appearing in heavy-ion collisions): several 10 -19 s. It is reasonable to introduce the 
evaporation prescription simply to the TDDFT final products. In this context the TDDFT final 
fragments have the meaning of products just after the early stage of heavy-ion reactions (several 
lO" 21 s). 
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Figure 2. (Colour online) Initial configuration for 248 Cm + 48 Ca; for instance, a case with 
b = 10 fm is depicted. Presented box is exactly the same as that used in calculations 
(48 x 48 x 24 fm 3 ). 



2.2. Evaporation prescription 

In complete fusion reactions the cross-section for the formation of a certain evaporation residue 
is usually given by three factors [5]: 

°ER{Ecm) =^2crCp(E cm , J) X P C N(E cm , J) X P S v(E cm , J) (1) 

J 

where ctcp, Pcn and Psv mean the capture cross-section, the probability for the compound 
nucleus formation, and the probability for survival of the compound nucleus against fission. All 
three factors are functions of the centre-of-mass energy E cm and the total angular momentum J, 
where J can be related with the impact parameter. For light systems Pcn and Psv are about 
unity and ger ~ J2j &CP- But in superheavy systems the strong Coulomb repulsion and large 
angular momenta lead to small values of Pcn and Psv and therefore to the small cross-sections of 
the evaporation residues observed in the experiments. That means, different from light systems, 
it is necessary to introduce additional thermal effects for the superheavy element synthesis. 
First, a cp is sufficiently considered in the TDDFT if we restrict ourselves to a sufficiently high 
energy exceeding the Coulomb barrier (cf. sub-barrier effects such as tunnelling are not taken 
into account in the TDDFT). Second, Pcn is fully considered in the TDDFT, which is a kind of 
mass equilibration also related to charge equilibration. Third, Psv whose relative time-scale is 
by no means equal to the former two probabilities is not satisfactorily considered in the TDDFT. 
This probability is much more related to thermal effects. Consequently, further consideration is 
necessary only for Psv as far as the energy above the Coulomb barrier is concerned. 

Several factors are included in Psv such as probabilities for fission of the compound nucleus, 
neutron-evaporation, proton-evaporation, deuteron-evaporation, alpha-particle-evaporation and 
so on. As a very first trial of our long-term project of having a generalization of TDDFT 
calculations, let us take the probability for neutron-evaporation and alpha-particle evaporation: 

PSV '•— (1 — Pn,evap)d ~ Pa,evap) 




Figure 3. (Colour online) Disappearance of charge equilibration for a higher incident energy 
(1692 MeV), where a case with b = 10 fm is shown at 18.6xl0 -22 s. The charge equilibration 
upper energy limit is calculated to be E/A — 3.7 MeV using the formula shown in Ref. [3], 
and this incident energy is E/A = 5.7 MeV. 19 (JV/Z=1.375) and the fissioning composite 
nucleus " 277 108" are obtained. Here the production of 19 whose N/Z is out of the initial range 
"1.40 < N/Z < 1.58" indicates the disappearance of charge equilibration. Note again that the 
periodic boundary condition is applied here. 

where this choice of neutron and alpha-particle evaporation is based on the experimental insight 
and results of the present reaction 248 Cm + 48 Ca. 

The algorithmic description of the evaporation prescription is as follows: 

1) obtain TDDFT fragments (Ai and Zj) with excitation (E*) and kinetic (Ki) energies; 

2) for each fragment calculate evaporation energy: 



where Ei^ evap and Ei^ nevap are the cooling energies due to one alpha-particle and one neutron 
emissions, respectively, and £^ ; ^ n and Ei^ n kin denotes the kinetic energy of alpha-particle and 
neutron, respectively (cf. Boltzmann distribution); 
3) for each fragment find the maximum integer n^ a satisfying 



where n^ a means the maximum number for alpha-particle emission from the i-th fragment; 
4) for each fragment find the maximum integer n^ n satisfying 



where n^ n means the maximum number for neutron emission from the z-th fragment. 

For the z-th fragment, a fragment with mass number Ai — 4n^ a — n^ n and Z{ — 2rii^ a is pro- 
duced by emitting n^ a times an alpha-particle and n^ n times a neutron. The excitation energy 
for the fragment is equal to E* — ni^ a Ei^ a ^ evap — ni^ n Ei^ n ^ evap . Consequently, the total number 
for alpha-particle and neutron emissions are determined by n a = an d n n — 




E* > n^ a Ei t 




Energy [MeV] 




| | Fusion 

| | Deep-inelastic 

□ Elastic 



4 6 8 10 12 14 

b[fm] 



Figure 4. (Colour online) Diagram of different reaction channels reactions obtained in 
248 Cm+ 48 Ca collisions by TDDFT. The preferred reaction channels for different beam energies 
and impact parameters are given. The beam energies are all located above the Coulomb barrier 
which is 209.00 MeV. The results show fusion, deep-inelastic, and elastic events. 



respectively. Hereafter we assume that neutrons, which are easier to be emitted compared to 
alpha-particles, are emitted earlier than alpha particles. This evaporation prescription is prefer- 
ably implemented directly to the TDDFT results without having the particle number projection 
or similar statistical treatment for the mass distribution of final fragments. Indeed, if such a 
treatment is applied, the excitation energy of the obtained fragments cannot be identified, and 
shell effects, which are expected to be included in the pure TDDFT results, might be spoiled. 



3. TDDFT calculations for superheavy element synthesis 

Time-dependent density functional calculations with a Skyrme interaction (SLy6 [6]) are carried 
out in a spatial box of 48 x 48 x 24 fm 3 with periodic boundary condition. The unit spatial 
spacing and the unit time spacing are fixed to 1.0 fm and 2/3 x 10 -24 s, respectively. 

The initial positions of 248 Cm and 48 Ca are fixed to (0,6,0) and (-15,0,0), respectively (Fig. 
[2]). The initial 248 Cm is almost spherical; the diameter for x, y, and z directions are 19 fm, 
19 fm and 18 fm, respectively. 248 Cm (the right hand side on x — y-plane) does not have initial 
velocity on the frame, while the initial velocity parallel to the x-axis is given to 48 Ca (the left 
hand side on x — y-plane) . 

First of all, the disappearance of charge equilibration for higher energies can be obtained 
in such reactions producing very heavy nuclei (Fig. [3]). This event is a fragmentation without 
having a compound nucleus, so that it is rather similar to quasi-fission event. In this paper, the 
difference between fusion-fission and quasi-fission is identified by whether the fission products 
satisfy "1.40 < N/Z < 1.58" or not. According to the disappearance of charge equilibration, the 
nucleus with its N/Z satisfying N/Z < 1.40 or N/Z > 1.58 can be obtained for higher energies, 
but the obtained fragments are light. For the synthesis of neutron-rich and heavy nuclei such as 
heavy elements located at the island of stability, there is a contradicting situation where high 
energy (disappearance of charge equilibration) is preferred for neutron-rich synthesis, while low 
energy (appearance of charge equilibration) is suitable for heavy synthesis. With respect to 
producing superheavy elements in the laboratory, here is a reason why we should not go further 
into higher energies. 

The systematic results of TDDFT calculations are summarized in Fig. 1 At a glance, 
these results, which include many fusion events, provide a quite optimistic view for producing 
superheavy elements. However, in comparison with experiments, the corresponding fusion cross- 
section of those low-energy heavy-ion reactions is too high to believe. Consequently, although 
these TDDFT results are still legitimate to show products just after the early stage of heavy- 




Figure 5. (Colour online) Time evolution of 248 Cm+ 48 Ca with a fixed incident energy 352 MeV 
in the centre-of-mass frame, where a box is fixed to 48x48x24 fm 3 . Snapshots at 2.3xl0 -22 s, 
8.7xl0" 22 s, 15.1xl0" 22 s and 21.5xl0" 22 s are shown, where a snapshot at 21.5x10 s is 
not shown only for b = 12 fm. Fusion, deep-inelastic, and elastic events appear depending on b. 
In particular, the life-time of composite nucleus for b = 10 fm is roughly equal to 15xl0 -22 s. 



ion reactions, it is necessary to take into account Psv m order to have comparable results to 
experiments. 

To compare the difference between the results before and after the evaporation prescription, 
let us take cases with the incident energy 352 MeV. TDDFT calculations with different b are 
shown in Fig. [5j Resulting reactions are 



248 Cm 


+ 


48 Ca - 


+ 296 Lv 




248 Cm 


+ 


48 Ca - 


+ 244 Cm+ 52 Ca 


(2) 


248 Cm 


+ 


48 Ca - 


+ 248 Cm+ 48 Ca 





for 6 = 8, 10, 12 fm, respectively. The first, second and third cases correspond to fusion, fission of 
the composite nucleus (charge equilibration appears), and elastic scattering, respectively. If the 
evaporation prescription is introduced to cases with the incident energy 352 MeV, the resulting 



reactions shown in Eq. ([2|) become 



248 Cm 


+ 


48 Ca -> 


275 Hs + 4a + 5n 




248 Cm 


+ 


48 Ca -> 


226 Th + 50 Ca + 3a + 8n 


(3) 


248 Cm 


+ 


48 Ca -> 


248 Cm + 48 Ca 





for 6 = 8, 10, 12 fm, respectively. Alpha-particles and neutrons are emitted for b < 10 fm. In the 
same manner, in case of the incident energy 436 MeV (the higher energy), the resulting reaction 
is 

248 Cm + 48 Ca ^ 264 Rf + g a + 8n 

for b < 10 fm. In case of the incident energy 268 MeV (the lower energy), the resulting reactions 
are 

248 Cm + 48 Ca ^ 280 Ds + 3a + 4n 

248 Cm + 48 Ca ^ 238 P U + 47 Ar + 2a + 3n 
248 Cm + 48 Ca ^ 237 Np + 48 K + 2a + 3n 
for b = 4, 6, 8 fm, respectively. 

4. Summary 

As a result, for low-energy reactions above the Coulomb barrier energy, the obtained reactions 
are summarized. 

E = 268MeV, 

248 Cm+ 48 Ca ^ 292 Lv + 4n 

-> 288 Fl + a + 4n 

-»■ 284 Cn + 2a + 4n 

-»■ 280 Ds + 3a + 4n, (6 < 4 fm) 

248 Cm+ 48 Ca ^ 246 Cf + 47 Ar + 3n 

-> 242 Cm + 47 Ar + a + 3n 

-»■ 238 Pu + 47 Ar + 2a + 3n, (6 = 6 fm) 

248 Cm+ 48 Ca ^ 2 45 Bk + 48 K + 3n 

-> 241 Am + 48 K + la + 3n 

-> 237 Np + 48 K + 2a + 3n, (6 = 8 fm) 

248 Cm + 48 Ca _^ 248 Cm + 48 Ca q > 1Q fm ). 



E = 352MeV, 

248 Cm+ 48 Ca ^ 291 Lv + 5n 

->• 287 Fl + a + 5n 

-»■ 283 Cn + 2o; + 5n 

-»■ 279 Ds + 3a + 5n 

-> 275 Hs + 4a + 5n, (6 < 8 fin), 

248 Cm+ 48 Ca ^ 238 Cm + 50 Ca + 8n 

-> 234 Pu + 50 Ca + a + 8n 
-»■ 230 U + 50 Ca + 2a + 8n 
-> 226 Th + 50 Ca + 3a + 8n, (6 = 10 fin), 

248 Cm + 48 Ca ^ 248 Cm + 48 Ca ( & > u fm ). 

£ = 436MeV, 

248 Cm+ 48 Ca ^ 288 Lv + 8n 

-»■ 284 Fl + a + 8n 

-»■ 280 Cn + 2a + 8n 

-> 276 Ds + 3a + 8n 

-> 272 Hs + 4a + 8n 

-»■ 268 Sg + 5a + 8n 

-»■ 264 Rf + 6a + 8n, (6 < 10 fin), 
248 Cm + 48 Ca 248 Cm + 48 Ca, (6 > 12 fm), 

where the indicated energy is the centre-of-mass energy and b is incremented by 2 fm. Transfer 
of a few nucleons contributing to charge equilibration can be noticed for E = 268 and 352 MeV 
cases. The experimental results reveal that fission of the compound nucleus occurs with very 
large probability, which will be included in future steps into the calculations. 
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